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' Abstract 

^ Functional limit theorem for continuous-time random walks (CTRW) are found 

in general case of dependent waiting times and jump sizes that are also position 
dependent. The limiting anomalous diffusion is described in terms of fractional 
dynamics. Probabilistic interpretation of generalized fractional evolution is given 
in terms of the random time change (subordination) by means of hitting times 
processes. 
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^ 1 Introduction 

Suppose (Xi,Ti), (X2,T2),... is a sequence of i.i.d. pairs of random variables such that 
Xi e Ti G R+ (jump sizes and waiting times between the jumps), the distribution of 
each {Xi, Ti) being given by a probability measure ip{dxdt) on R'' x R+. Let 



Nt = max{n : T, < t}. 



i=l 



The process 

S^^=X, + X2 + ... + X^^ (1) 

is called the continuous time random walk (CTRW) arising from i/j. These CTRW were 
introduced in [T7] and found numerous applications in physics and economics (see e.g. 
[23] . [H], [3], [12], [15] and references therein). Of particular interest are the situations. 
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where Tj belong to the domain of attraction of a /? G (0, l)-stable law and Xj belong to the 
domain of attraction of a a e (0, 2) -stable law. The limit distributions of appropriately 
normalized sums were first studied in [7] in case of independent Tj and Xj (see also 
[TIJ). In [5] the rate of convergence in double array schemes was analyzed and in [Hj 
the corresponding functional limit was obtained, which was shown to be specified by a 
fractional differential equations. The importance of the analysis of the case of dependent 
Tj and Xj was stressed both in and [H]. Here we address this issue. Moreover we 
extent the theory to include possible dependence of (T„, X„) on the current position. Our 
method is quite different from those used in [7j, |TT], [M]. It is based on the finite difference 
approximations to continuous-time operator semigroups and applies the previous results 
of the author from |8j on stable-like processes. 

It was noted in [H] that fractional evolution appears from the subordination of Levy 
processes by the hitting times of stable Levy subordinators. Implicitly this idea was 
present already in [19]. We are going to develop here the general theory of subordination of 
Markov processes by the hitting time process showing that this procedure leads naturally 
to (generalized) fractional evolutions. In particular, in spite of the remark from [Hj that 
the method from (going actually back to [T7j) "does not identify the limit process" we 
shall give a rigorous probabilistic interpretation of the intuitively appealing (but rather 
formal) calculations from [19]. 

In the next Section we demonstrate our approach to the limits of CTRW by obtaining 
simple (but nevertheless seemingly new) limit theorems for position depending random 
walks with jump sizes from the domain of attraction of stable laws. In Section 3 these 
results will be extended to double scaled random walks, which are needed for the analysis 
of CTRW. Section 4 (which is independent of Section 2 and seems to be of independent 
interest) is devoted to the theory of subordination by hitting times. In Section 5 we 
combine the two bits of the theory from Sections 3 and 4 giving our main results on 
CTRW. 

Let us fix some (rather standard) notations to be used throughout the paper. For a 
locally compact space X we denote by C(X) the Banach space of bounded continuous 
functions (equipped with the the sup-norm) and by Coo(X) its closed subspace consisting 
of functions vanishing at infinity. We denote by (/, /i) the usual pairing J f{x)fi{dx) be- 
tween functions and measures. By a continuous family of transition probabilities (CFTP) 
in X we mean as usual a family p{x; dy) of probability measures on X depending con- 
tinuously on X G X, where probability measures are considered in their weak topology 
{fin — > /i as n — > oo means that (/, /x„) — >■ (/, /i) as n ^ oo for any / G C(X)). 

For a measure fi{dy) in R'^ and a positive number h we denote by fi{dy/h) the scaled 
measure defined via its action 



on functions g G C(R'^). 

The capital letters E and P are reserved to denote expectation and probability. The 
function 6{x) is the usual Dirac function (distribution). 
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2 Limit theorems for position dependent random walks 



For a vector y G we shall always denote by y its normalization y = y/\y\, where \y\ 
means the usual Euclidean norm. 

Fix an arbitrary a G (0,2). Let S : R'' x S'^^^ ^ R+ be a continuous non-negative 
function that is symmetric with respect to the second variable, i.e. S{x,y) = S{x, —y). 
It defines a family of a-stable (i-dimensional symmetric random vectors (depending on 
X G R'') specified by its characteristic function (p^ with 



J s-^ 



In 0.(p) = / / . ( e^^-'^^ - 1 - J ^^(^' dsi. (2) 



where ds denotes the Lebesgue measure on the sphere S*^ ^ . It is well known that it can 
be also rewritten in the form 

ln0.(p)=C„ / \{v.lTS{x,i)dsi 
with a certain constant Ca- 

Remark 1 There are no obstacles for extending our theory to non-symmetric stable laws. 
But working with symmetric laws shorten the formulas essentially. 

Theorem 2.1 Assume 



Ci< / \iP,sTS{x,s)dss<C2 

for all p with some constants Ci, C2 and that S{x, s) has bounded derivatives with respect 
to X up to and inclusive order q > 3 (if a < 1, the assumption q > 2 is sufficient). Then 
the pseudo-differential operator 

Lf{x)=\nM-^)fix)= r [ {f{x + y)-f{x))-^S{x,y)dsy (3) 

^ ex Jo J s'^-^ \y\ 

generates a Feller semigroup Tt in Coo(R'^) with the space C^^^^BJ^) fl Coo(R'^) being its 
invariant core. 

This result is proven in [8] and [9]. 

Remark 2 In it is also shown that this semigroup has a continuous transition density 
(heat kernel), but we do not need it. 

Denote by Z.j.{t) the Feller process corresponding to the semigroup Tf. We are inter- 
ested here in discrete approximations to Tt and Zr^it). 
We shall start with the following technical result. 

Proposition 2.1 Assume that p{x]dy) is a CFTP in R"^ from the normal domain of 
attraction of the stable law specified by ([2]). More precisely assume that for an arbitrary 
open Q G 3"^^^ with a boundary of Lebesgue measure zero 

I I P{x]dy)r^-^f S{x,s)dss, n 00, (4) 

J\y\>n Jyen Ju 
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(i.e. the ratio of the two sides of this formula tends to one as n oo) uniformly in x. 
Assume also that p{x, {0}) = for all x. Then 

mm(l, \y\^)p{x, dy/h)h~" min(l, \y\'^)-^j^S{x, y)dsy, h ^ 0, (5) 

where both sides are finite measures on \ {0} and the convergence is in the weak sense 
and is uniform in x & R'^. If a < 1, then also 

71 I p 

mm{l, \y\)p{x, dy /h)h~'' ^mm{l,\y\)-—— S{x,y)dsy, h 0, 

\y\ Jn 

holds in the same sense. 

Remark 3 As the limiting measure has a density with respect to Lebesgue measure, the 
uniform weak convergence means simply that the measures of any open or closed set con- 
verge uniformly in x. 



Proof. By 



p{x;dz/h)h " = [ [ p{x;dy)h "~ — — f S{x,s)dss 

z\>AJzeQ J\y\>A/hJy&n ChA'^Jq 



as /i ^ 0. Hence 



p{x;dz/h)h / / S{x,s)dss. 



Rence p{x;dz/h)h ° converges weakly to \z\ ^°'~^^^d\z\S{x,z/\z\)ds{z/\z\) on any set sep- 
arated from the origin. It is easy to see that follows now from the uniform bound 

[ mm{l, \y\^)p{x, dy / h)h-'' KCe^-"" (6) 

with a constant C. In order to prove ([6]) let us observe that 

/ p{x, dy) < Cn-" 

J\y\>n 

with a constant C uniformly for all x and n > (in fact it holds for large enough n by 
IHand is extended to all n, because all p{x,dy) are probability measures). Hence for an 
arbitrary e < 1 one has 

[ mm{l, \y\'^)p{x, dy / h)h-'' = I h'^\z\'^)p{x,dy/h)h-''. 

J\y\<e J\z\<e/h 

Representing this integral as the countable sum of the integrals over the regions 

e/{2''+^h) <y< e/{2^h), 

it can be estimated by 

oo 2 °° 

k=0 k=0 



This yields ([6]), since the sum on the r.h.s. converges. 

The improvement concerning the case a < 1 is obtained similarly. 
Consider the jump-type Markov process Z^(t) generated by 

(L,/)(x) = {fix + hy) - f{x))p{x- dy) (7) 

For each h the operator is bounded in Coo(R'^) and hence specifies a Feller semigroup 
there. The probabilistic interpretation of Z^{t) is as follows. Starting at a point x one 
waits a random 6 = /^""-exponential time r (i.e. distributed according to P(r > t) = 
exp{—t6)) and then jumps to x -|- hY, where Y is distributed according to p{x] dy). Then 
the same repeats starting from x + hY, etc. In case when p does not depend on x 

Z\t) = h{Yi + ... + YNj 

is a normalized random walk with the number of jumps A^^ being a Poisson process with 
parameter so that ENt = th~°'. In particular, the number of jumps n = Nt ^ th'"" 
for small h so that Z^{1) ~ + ... + Y^). 

Theorem 2.2 The semigroup T/* generated hy converges to the semigroup Tt generated 
hy L. In particular, the corresponding processes converge in the sense of finite- dimensional 
marginal distrihutions. 

Remark 4 Everywhere in this paper we work with the convergence of semigroups only. 
However hy the standard results (see e.g. Theorem 19.25 in 16]) for Feller processes 
this convergence is equivalent to the convergence of the distrihutions of trajectories in an 
appropriate Skorokhod space of cddldg paths. 

Proof By (JTD 

= ^J {fix + z)~ f{x))p{x; dz/h), 

and by Proposition 12. II this converges to Lf[x) as /i — > uniformly in x for / G Coo(R'^) H 
C^(R'^). By a well known result (see e.g. [I3]) the convergence of the generators on the 
core of the limiting semigroup implies the convergence of semigroups. 

The next result concerns the approximations with a non-random number of jumps. 
Define the process S^{t) = 5'^([t]) (by the square bracket the integer part of a real number 
was denoted) via 

S^{0) = x, S^^il) = X + hY,, S',{j)=S'Aj-^) + hY„... 

where each Yj is distributed according to p{Sj-i,dy). If p{x;dy) does not depend on x, 
then 

Sf^{n)=x + hiY, + ... + Yn) 

is just a standard random walk. 

We like to compare the Feller process Zx{t) on an arbitrary fixed time interval [0,to] 
with the discrete approximations S^{t/T), when the number of jumps n = t/r is connected 
with the scaling parameter h hj t = h°'. 
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Theorem 2.3 Under the assumptions of Theorem \2.1\ and Proposition \2. 1\ for any f G 
Coo(R'^), Ef{Sl^{t/T)) converges to Ttf[x) uniformly on t & [0,to]; as r = /i" — > 0. 
In particular, the processes Sl^{t/T) converge to Zx{t) in the sense of finite- dimensional 
distributions. 

Proof. It is enough to prove the required convergence for / G C^(R'^) fl Coo(R'^) only (by 
Theorem 12. ip . Let such an / be chosen. Denote fk{x) = Ef{S!^{k)). Then by the Markov 
property fk = Rhf, where the operator is defined via the formula 

Rhf{x) = j f{x + hy)p{x; dy). 

Clearly each is a positivity preserving contraction on Coo(R"'). On the other hand, 
the recurrent equation = Rfk-i can be rewritten as 

fk{x)-f,.,{x) ^ r ^ _ (8) 



And this is a discrete time approximation to the equation 

on the functions / G C^(R'^) fl Coo(R'^) (and differentiable in t). Since this scheme is 
well-posed and stable (as it is solvable uniquely by the contraction i?^) and the solution 
to ([9]) is uniquely defined and preserves the space C^(R'^) fl Coo(R'^) (by Theorem 12.11) . 
it follows by the standard (and easy to prove) general results (see e.g. [20]) that the 
solutions to the finite-difference approximation converge to the solution of ([9]). Theorem 
is proved. 

In case of p not depending on x, Theorem 12.31 turns to the known fact on the con- 
vergence of random walks with the distribution of jumps from the domain of normal 
attraction of a stable law to the corresponding stable Levy motion. 

3 Double-scaled random walks 

To apply the developed theory to CTRW we shall need a generalization with multi-scaled 
walks that we present now. 

We are interested in a process in R*^ x R_|_ specified by the generator 

Cf{x,u)= / {f{x + y,u) - f{x,u)) S{x,u,y)dsy 
Jo Js'i-^ \y\~^ 

f°° 1 
+ J {f{x,u + v) - f{x,u))-^w{x,u)dv. (10) 

The following result (and its proof) is a straightforward generalization of Theorem 12. 1[ 

Theorem 3.1 Assume 

Ci< j \{p,s)\''S{x,u,s)dss<C2, Ci<w{x,u)<C2 
Js'i-^ 
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with some constants Ci,C2 and that S{x,s) and w{x,u) have bounded derivatives with 
respect to x and u up to and inclusive order q > 3. Then the pseudo- differential operator 
(fTOj) generates a Feller semigroup % in Coo(R'^ x R+) (continuous functions up to the 
boundary) with the space {C^^^ fl Coo)(R'^ x R_|_) being its invariant core and hence a 
Feller process {Y, V) (t) in BJ^ x R_|_ . 

We shall obtain now the corresponding extension of Theorems 12.2^ 12.31 

Theorem 3.2 Assume p{x,u; dydv) is a CFTP in R*^ x R_,_, which is symmetric with 
respect to the reflection y i-^ —y and for which 

p{x, u] {0} X R+) + p{x, u] R'^ X {0}) = 0. 

Assume also that the projections belong to the domain of normal attraction of stable laws; 
more precisely, that uniformly in {x, u) 

I I p{x,u;dydv) — - f S{x,u,s)dss, n ^ oo, (11) 

J\y\>n Jyen Q^^" Ju 

and 

/ p{x,u; dydv) ——gW{x,u, A), n oo, (12) 

v>n J\y\>A 

for any A >0 with a measurable function w of three arguments such that 

w{x,u,0) = w{x,u). Mm w{x,u. A) = (13) 

A^oo 

(so that w{x, u. A) is a measure on R+ for any x,u). 
Consider the jump-type processes generated by 

(Crf) {x, u) = - [ {fix + r'/^y, u + r^'^v)) - f{x, u))p{x, u; dydv). (14) 



Then the Feller semigroups in Coo(R'^ x R_,_) of these processes (which are Feller, 
because Ch is bounded in Coo(R'^ x R+) for any h) converge to the semigroup %. 

Proof. As in Proposition 12. II one deduces from (fTTj) . (fT2|) that uniformly in x,u 

rtQo ^ I I 

min(l,|?/n/ p{x, u]dy/hdv)h-" ^ mm{l, \y\'^)-——S{x,y)dsy, h ^ 0, (15) 

Jo \y\ 

and 

f dv 

min(l,t>) / p{x,u; dydv/h)h^^ ^ mm{l,v)w{x,u, A)^—, /i — > 0, (16) 

J\y\>A 

Next, assuming / G (C^ fl Coo(R'^ x R_|_) and writing 

Crf{x,u) =1 + 11 

with 

/ = ^ / {f{x+T^^"'y,u)-f{x,u))p{x,u;dydv) + ^ j {f{x,u+T^"^v)-f{x,u))p{x,u;dydv) 



and 

//= [{f{xW^''y,u+T^/'^v)-f{xW^"y,u))-{f{x,uW^'^v)-f{x,u))]p{x,u;dydv) 

one observes that, as in the proof of Theorem \2.2\ f|T5l) and f|T6l) (the latter with A = 0) 
imply that / converges to Cf{x, u) uniformly in x, u. Thus in order to complete our proof 
we have to show that the function // converges to zero, as r — *• 0. We have 

II = J {g{x + r^^^y,u,v) — g{x^u^v))p{x^u\dydv /t^^^) — 

with 

g{x, M, v) = f{x, u + v)- f{x, u). 

By our assumptions on / 

df 

\g{x, u,v)\ < C min(l, t>)(max | — | + max |/|) < C min(l, v), 

and 

dg d'^f df 

It— (x, M, f )| <C min(l, t>)(max L „ I + max I7— 1) < C min(l, v) 
ox ouox ox 

with some constants C and C . Hence by f|T6l) and f|T3l) for an arbitrary e > there exists 
a A such that 

{g{x + T^''^y,u,v) - g{x,u,v))p{x,u]dydv/T^^^)- < e; 

y\>A T' 

and on the other hand, for an arbitrary A 

{g{x + T^^^y, u, v) — g{x, u, v))p{x, u; dydv/r^^^)- < t^^"Ak 

y\<A ^ 

with a constant k so that II can be made arbitrary small by first choosing large enough 
A and then choosing small enough r. 

Define now the process {Y, = (^) ^)^,m([V'^])> where 

{y.V)lu{^) = (F,V^)^,„(l) = (x + ri/-Fi,« + ri/^rO,..., 

{y. V)IM = (y, V)Uj - 1) + (ri/"F„ r'/f'V,), ... 

and each pair {Yj, Vj) is distributed according to p{(Y, V)l „(j — 1); dydv). If p(x, u; dydv) 
does not depend on x, m, then 

{y^ y)lM = {x,u) + (ri/"(yi + ... + r„),ri/^(\/i + ... + K)). 

In view of Theorem 13. 2l the following result is obtained by literally the same arguments 
as Theorem 12. 3[ 

Theorem 3.3 Under the assumptions of Theorems \ 3.1\ and \3.S\ the linear contractions 
Ef{{Y,V)l.^(t/T)) in Coo(R'^ x R^.) converge to the semigroup Ttf{x,u) of the process 
(Y,y)(t) uniformly ont & [0,to]; as r — 0. 
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4 Subordination by hitting times and generalized frac- 
tional evolutions 

Let X{u), M > be a Levy subordinator, i.e. an increasing i.i.d. cadlag Feller process 
(adapted to a filtration on a suitable probability space) with the generator 

Af{x) = l^^ifix + y)- f{x)Mdy) + a^, (17) 

where a > and u is a Borel measure on {y > 0} such that 

m.m{l, y)h'{dy) < oo. 



We are interested in the inverse function process or the first hitting time process Z{t) 
defined as 

Zxit) = Z{t) = mi{u : X{u) > t} = snp{u : X{u) < t}, (18) 

which is of course also an increasing cadlad process. To make our further analysis more 
transparent (avoiding heavy technicalities of the most general case) we shall assume that 
there exist e > and j3 G (0, 1) such that 

i^{dy)>y'+^, 0<y<e. (19) 

For convenient reference we collect in the next statement (without proofs) the elementary 
(well known) properties of X{u). 

Proposition 4.1 Under condition ( fT9l) (i) the process X{u) is a.s. increasing at each 
point, i.e. it is not a constant on any finite time interval; (ii) distribution of X{u) for 
u > has a density G{u,y) vanishing for y < 0, which is infinitely differentiable in both 
variable and satisfies the equation 

^^A'G, (20) 

where A* is the dual operator to A given by 

A* fix) = j^U\^ -y)- f{x)Hdy) - a^, 

(Hi) if extended by zero to the half-space {t < 0} the locally integrable function G{t,y) on 
specifies a generalized function satisfying (in the sense of distribution) the equation 

— = A^G^b{u)b{y). (21) 

Corollary 1 Under condition f|T9l) (i) the process Z{t) is a.s. continuous and Z{0) = 0; 
(ii) the distribution of Z{t) has a continuously differentiable probability density function 
Q{t,u) for u > given by 

Q{t,u) = -^J^ G{u,y)dy. (22) 
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Proof, (i) follows from Proposition 14.11 (i) and for (ii) one observes that 

POO pt 

P{Z{t) <u)= P{X{u) >t)= Giu, y)dy=l - / Giu, y) dy 

Jt Jo 

which imphes ( 122|) by the different iabihty of G. 



Theorem 4.1 Under condition (1191) the density Q satisfies the equation 

A^Q = ^ (23) 

for u > 0, where A* acts on the variable t, and the boundary condition 

\im Q{t, u) = -A* e{t) (24) 

where 9{t) is the indicator function equal one (respectively 0^ for positive (respectively 
negative) t. If Q is extended by zero to the half-space {u < 0}, it satisfies the equation 

A*Q = ^ + 5{u)A^e{t), (25) 
ou 

in the sense of distribution (generalized functions). 

Moreover the (point-wise) derivative ^ also satisfies equation (!23l) for u > and 
satisfies the equation 

in the sense of distributions. 

Remark 5 In the case of a j3-stable subordinator X{u) with the generator 

^/W = -Fr-^/ U{^ + y)-f{^))y-''^dy, (27) 



r(-/9) Jo 

one has 

(these equations can be considered as the definitions of fractional derivatives; we refer to 
books IW^ and IW^ for a general background in fractional calculus; a short handy account 
is given also in Appendix to JT^]), in which case equation fl25p takes the form 

d^Q dQ , t-^ 

d^ + ^ = ^(")f(r^ ^''^ 

coinciding with (BI4) from [T^ . 

Proof. Notice that by fl2^ . fl2Ul) and by the commutativity of the integration and A* one 
has 

pi o pt, pt 

Q{t,u) = - Giu,y)dy=- {A^G{u, .)){y) dy = -A^ G{u,y)dy. 



00 du 
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This implies ( 123|) (by differentiating with respect to u and again using (!22l) ) and fl24l) . 
because G(0,?/) = 

Assume now that Q is extended by zero to {m < 0}. Let be a test function (infinitely 
differentiable with a compact support) in R^. Then in the sense of distribution 

;|;-A*w>) = (o.(-|;-^w 

/■oo /■ g 

= lim / du dtQ(t,u)(-- A)(j)(t,u) 

Je Jn du 

= lim f du f dt(f)(t,u)(-^ - A*)Q(t,u) + [ (j)(t,e)Q(t,e) dt 
Ue Jn du 

The first term here vanishes by ( l23l) . Hence by (E 



which clearly implies fl25l) . The required properties of ^ follows similarly from the 
representation 

dQ dG 
^(t,«) = -^(n,t). 

We are interested now in the random time change of Markov processes specified by 
the process Z{t). 

Theorem 4.2 Under the conditions of Theorem 4-1 ^^t Y{t) be a Feller process in R"^, 
independent of Z(t), and with the domain of the generator L containing {Coo nC^)(R'^). 
Denote the transition probabilities ofY{t) by 

T{t,x,dy) = P{Y,it) G dy) = P.(F(t)). 

Then the distributions of the (time changed or subordinated) process Y{Z{t)) fort > 
are given by 

poo 

P^{Y{Z{t)) e dy) = T{u,x,dy)Q{t,u)du, (30) 
Jo 

the averages f(t,x) = Ef(Yx{Z{t))) of f E {Coo H C^)(R"') satisfy the (generalized) 
fractional evolution equation 

Atf{t, x) = -LJ{t, x) + f{x)A*e{t) (31) 

(where the subscripts indicate the variables, on which the operators act), and their time 
derivatives h = df/dt satisfy for t > the equation 

A*h = -L,h + f{x)j^A*9{t). (32) 

Moreover, ifY{t) has a smooth transition probability density so thatT{t, x, dy) = T{t, x, y)dy 
and the forward and backward equations 



dT 
'dt 



{t, X, y) = L^T{t, X, y) = LlT{t, x, y) (33) 
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hold, then the distributions ofY{Z(t)) have smooth density 

POO 

9{t,x,y)= T{u,x,y)Q{t,u)du (34) 

satisfying the forward (generalized) fractional evolution equation 

A*g = ~Llg + 6{x - y)A*9it) (35) 
and the backward (generalized) fractional evolution equation 

A^g = -L,g + 6{x - y)A*e{t) (36) 
with the time derivative h = dg/dt satisfying for t > the equation 

A*h = -L*h + 5{x-y)^A*e{t) (37) 
^ at 

Remark 6 In the case of a jS -stable Levy subordinator X{u) with the generator fl27|) . 
where fl28l) hold, the left hand sides of the above equations become fractional derivatives 
per se. In particular, if Y(t) is a symmetric a-stable Levy motion, equation ( l35l) takes 
the form 

-— To(t, y — x) = ^, , q(t, y — x) + S(y — x)— --, (38) 

deduced in fW^ and l21^ . The corresponding particular case of also appears in JT^F 
as well as in JJE/, where it is called a formula of separation of variables. Our general 
approach makes it clear that this separation of variables comes from the independence of 
Y{t) and the subordinator X{u) (see Proposition 4^ for a more general situation). 

Proof. For a continuous bounded function / one has for t > that 

poo POO 

Ef{Y^{Z{t)) = / E{f{Y^{Z{t))\Z{t)=u)Q{t,u)du= / Ef{Y^{u))Q{t,u)du 
Jo Jo 



by the independence of Z and Y. This imphes fl30|) and flMl) . 
From Theorem 14.11 it follows that for t > 

/CO POO Q 

G{u, x,y)A*Q{t,u) du = lim J G{u, x,y)—Q{t,u) du 

roo Q 

= - J -Q-G{u, X, y)Q{t, u) du + 5{x - y)A*9{t), 

where by fl33l) the first term equals —L^g = L^g, implying fl35l) and fl36|) . Other equations 
are proved analogously. 

Now we like to generalize this theory to the case of Levy type subordinators X{u) 
specified by the generators of the form 

POO gn 

^/(a;) = y^ ifix + y)-fix)Mx,dy) + aix)^ (39) 

with position depending Levy measure and drift. We need some regularity assumptions 
in order to have a smooth transition probability density like in case of the Levy motions. 
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Proposition 4.2 Assume that (i) v has a density v{x, y) with respect to Lebesgue measure 
such that 

Cimin {y'^''^\y~^~^') < iy{x,y) < ^2 max {y'^'^\y'^~^') (40) 

with some constants Ci,C2 > and < (3i < (32 < 1 (ii) v is thrice continuously dif- 
ferentiable with respect to x with the derivatives satisfying the same estimate ( HOjl . (Hi) 
a{x) is non-negative with bounded derivatives up to the order three. Then the generator 
( !39l) specifies an increasing Feller process having for u > a transition probability den- 
sity G{u,y) = P{X{u) e dy) (we assume that X{u) starts at the origin) that is twice 
continuously differentiable in u. 

Remark 7 Condition (j40|) holds for popular stable-like processes with a position depen- 
dent stability index. 

Proof. The existence of the Feller process is proved under much more general assumptions 
in [Tj. A proof of the existence of a smooth transition probability density is given in [8] 
under slightly different assumptions (symmetric multidimensional stable-like processes), 
but is easily seen to be valid in the present situation. 

One can see now that the hitting time process defined by (ITSl) with X{u) from the 
previous Proposition is again continuous and has a continuously differentiable density 
Q(t,u) for t > given by (122]) . However (1251) does not hold, because the operators A 
and integration do not commute. On the other hand, equation (l26l) remains true (as 
easily seen from the proof). This leads directly to the following partial generalization of 
Theorem 14. 2[ 

Proposition 4.3 Let Y{t) be the same Feller process in R'^ as in Theorem \4.3\ but in- 
dependent hitting time process Z{t) be constructed from X{u) under the assumptions of 
Proposition \4.2\ 

Then the distributions of the (time changed or subordinated) process Y{Z{t)) for 
t > are given by (l30l) and the time derivatives h = df/dt of the averages f{t,x) = 
Ef(Yx{Z(t))) of continuous bounded functions f satisfy (ITTI) . 

At last we like to extend this to the case of dependent hitting times. 

Proposition 4.4 Let {Y,V){t) be a random process in R'^ x R+ such that (i) the com- 
ponents Y{t),V{s) at different times have a joint probability density 

(j){s, U] y, v) = P{Y{s) G dy, V{u) G dv) 

that is continuously differentiable in u for u,s > 0, and (ii) the component V{t) is in- 
creasing and is a.s. not a constant on any finite interval. For instance, the process from 
Theorem \3.1\ enjoys these properties. Then (i) the hitting time process Z{t) = Zyif) (de- 
fined by (|T8l) with V instead of X) is a.s. continuous, (ii) there exists a continuous joint 
probability density ofY{s),Z(t) given by 

d f°° 

gYis),zit){y,u) = — J (f){s,u;y,v)dv (41) 
and (Hi) the distribution of the composition Y{Z{t)) has the probability density 

/oo poo / Q poo \ 

gY{s),z{t){y,s)ds = J (f){s,u;y,v)dvj\u=sds. (42) 

Proof, (i) and (ii) are straightforward extensions of the Corollary to Proposition 14. 1[ 
Statement (iii) follows from conditioning and the definition of the joint distribution. 



13 



5 Limit theorems for position dependent CTRW 

Now everything is ready for our main result. 

Theorem 5.1 Under the assumptions of Theorems \3.1\ and \3.2\ let Z'^{t), Z{t) be the 
hitting time processes for V^lt/r) and V{t) respectively (defined by the corresponding 
formula flTHl) ). Then the subordinated processes Y'^{Z'^{t)/T) converge to the subordinated 
process Y{Z{t)) in the sense of marginal distributions, i.e. 

E,^of{Y^{Z^it)/T)) ^ E^^^iYiZm, r ^ 0, (43) 

for arbitrary x G R'^, / G Coo(R'^ x R+), uniformly for t from any compact interval. 

Remark 8 We show the convergence in the weakest possible sense. It does not seem dif- 
ficult to extend it to the convergence in the Skorokhod space of trajectories using standard 
tools (compactness etc) or the theory of continuous compositions from [2^]. Similar result 
holds for the continuous time approximation from Theorem \3.2[ 

Proof. Since the time is effectively discrete in V^lt/r), it follows that 

Z^{t) = max{u : X{u) < t}, 

and that the events {Z'^{t) < u) and {y'^{u/T) > t) coincide, which implies that the con- 
vergence of finite dimensional distributions of (F'^(s/r), V^i^u/r)) to (Y{s), V{u)) (proved 
in Theorem 13.31) is equivalent to the corresponding convergence of the distributions of 
(y-(s/r),Z-(t)) to (F(s),Z(t)). 

Next, since V^(0) = 0, is continuous and V{u) — * oo as m ^ oo and because the 
limiting distribution is absolutely continuous, to show (H3l) it is sufficient to show that 

P^,o[Y^{Z],{t)/r) eA]^ P^,o[Y{ZK{t)) G A], r ^ 0, (44) 

for large enough K > and any compact set A, whose boundary has Lebesgue measure 
zero, where 

ZJ^{t) = Z^t), < Z^t) < K, 

and vanishes otherwise, and similarly Zxit) is defined. 
Now 

K/t 

P[Y^{ZJ,{t)/T)eA]= J2 P[V^k)eAkZ^{t)e[kT,{k + l)T)] (45) 

k = l/KT 

and 

P[Y{ZK{t))eA]= gYis),zit){y,s)ds, (46) 

k=l/ Kt 

which can be rewritten as 
Y j dy j 9Y(Tk),z(t){y,s)ds+ ^ dy {gY(s),z(t)-9Y{Tk),z{t)){y, s) ds. 

k=l/KT-^^ •^'^^ k=l/KT-^^ •''^^ 

(47) 
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The second term here tends to zero as r — > due to the continuity of the function (HTI) . and 
the difference between the first term and (H5|) tends to zero, because the distributions of 
(F^(s/r), Z''{t) converge to the distribution of {Y{s), Z{t)). Hence (jHj) follows. Theorem 
is proved. 

In the case when S does not depend on u and w does not depend on x in (ITUl) . 
the limiting process {Y, V) (t) has independent components so that the averages of the 
limiting subordinated process satisfy the generalized fractional evolution equation from 
Proposition l4.3l and if moreover w is a constant, they satisfy the fractional equations from 
Theorem 14.21 In particular, if p{x, u, dydv) does not depend on (x, u) and decomposes 
into a product p{dy)q{dv), and the limit V{t) is stable, we recover the main result from 
[Hj (in a slightly less general setting, since we worked with symmetric stable laws and 
not with operator stable motions as in ^MjI), as well as of course the corresponding results 
from [7], [11] (put t = 1 in (H3!) ) on the long time behavior of the normalized subordinated 
sums ([1]). 
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